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ABSTRACT 

Highly ionized species such as C IV, N V, and O VI, are commonly observed in diffuse gas in various 
places in the universe, such as in our Galaxy's disk and halo, high velocity clouds (HVCs), external 
galaxies, and the intergalactic medium. These ions are often used to trace hot gas whose temperature 
is a few times 10^ K. One possible mechanism for producing high ions is turbulent mixing of cool gas 
(such as that in a high or intermediate velocity cloud) with hotter (a few times 10^ K) gas in locations 
where these gases slide past each other. By using hydrodynamic simulations with radiative cooling 
and non-equilibrium ionization (NEI) calculations, we investigate the physical properties of turbulent 
mixing layers and the production of high ions (C IV, N V, and O VI). We find that most of the mixing 
occurs on the hot side of the hot/cool interface, where denser cool gas is entrained and mixed into 
the hotter, more diffuse gas. Our simulations reveal that the mixed region separates into a tepid zone 
containing radiatively cooled, C IV-rich gas and a hotter zone which is rich in C IV, N V, and O VI. 
The hotter zone contains a mixture of low and intermediate ions contributed by the cool gas and 
intermediate and high-stage ions contributed by the hot gas. Mixing occurs faster than ionization or 
recombination, making the mixed gas a better source of C IV, N V, and O VI in our NEI simulations 
than in our collisional ionization equilibrium (CIE) simulations. In addition, the gas radiatively cools 
faster than the ions recombine, which also allows large numbers of C IV, N V, and O VI ions to linger 
in the NEI simulations. For these reasons, our NEI calculations predict more C IV, N V, and O VI 
than our CIE calculations predict. We also simulate various initial configurations and find that more 
C IV is produced when the shear speed is smaller or the hot gas has a higher temperature. We find 
no significant differences between simulations having different perturbation amplitudes in the initial 
boundary between the hot and cool gas. We discuss the results of our simulations, compare with 
observations of the Galactic halo and highly ionized HVCs, and compare with other models, including 
other turbulent mixing calculations. The ratios of C IV to N V and N V to O VI are in reasonable 
agreement with the averages calculated from observations of the halo. There is a great deal of variation 
from sightline to sightline and with time in our simulations. Such spatial and temporal variation may 
explain some of the variation seen among observations. 

Subject headings: Galaxy: halo — hydrodynamics — methods: numerical — turbulence — ultraviolet: 
ISM 



1. INTRODUCTION 

Highly ionized species such as C IV, N V, and 
O VI are often used as tracers for diffuse gas within 
the temperature range of (1 — 3) x 10^ K. These 
high ions are found by absorption line measurements 
in various pla.ces i n the universe, inclu d ing the disk 
(iBowen et all l2008t ISallmen et all 120081: iSavaee et al.l 



2001at 



JenkinsI 



Savage fc Massal 



1978 



an d 
'2005'; 



I1987t ICowieetal.1 I1981t 
the halo of the Milky Way 
'Sava ge et al.l 120031: iZsareo et al.l 



jGangulv et al 

2003; Ste rling et al _ 

1997; Se mbach et al.l ll997nSavage fc SembachI 11994 : 
Sembach fc Savage! Il9"m ~ The emission lines of these 



2 002t ISavage et all 12'ooTb 



high i ons ar e also found i n the Galaxy (iShelton et al.l 
20011 [20071 I Dixon et all [20M \Otte fc DixonI 1200ft 
Korpela et ain 2006: Welsh et al.ll2007| ). 
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and halo 
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observed high 
nearly stationary. 

There is also a population of high ions moving at veloc- 
ities of a few hundreds of k m s~^ with respect to the lo- 
cal standard of rest CLSR) (iCoUins et al.l l2007l: iFox et al.l 
[2005l[200llSembach et al1l2003n . The distances to some 
of these ions are unknown. Most of these ions are on 



sightlines that intersect H I high velocity clouds (HVCs) , 
but some are not. Thus, these latter ions are due to 
ionized HVCs. 

High ions are also found in exter nal galaxies via 
absorption line measurement s (LMC: [ Lehner fc Howld 
[20071: IDanforth fc BlaE [20Q1 ISankrit et al.l I2004D and 
O VI emissio n measurements (fBregman et al.l l2006bl : 
iGangulv et all 12006). High ions are detected not only 
in nearby galaxies b ut also in damped Lym an- a systems 
for distant galaxies ([Fox et al.l 120091 12007D and gamma- 
ray burst (GR B) host galaxies ([Prochaska et al.l 120081 : 
iFox et al]|2008f) . Even the intergalactic medium contains 
high ions. In the nearby intergalactic space, high ions are 
detec ted along sightlines to distant QSO's ([Tripp et al.l 
l2008f ). In the cooling flows of clusters of galaxies such 
as Abeh 426, Abeh 1795, and Abell 2597, high ions 
are observed in the flow gas cooling from hotter. X-ray 
emitting gas ([Bregman et al.ll2006ai lOegerle et al.|[2001l : 
IDixon et al.lll9960 . 

Observations of G IV, N V, and O VI in various places 
in the universe indicate that these ions are very com- 
monly produced. How they are produced and what their 
existence implies about the local physical conditions have 
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been longstanding questions. For example, in order to 
trace the detailed physical condition of the plasma in 
the Galactic disk and halo, ratios between the quanti- 
ties of different ions, either their absorption column den- 
sity or their emission intensity, are measured and com- 
pared with various model p redictions (Gnat & Ste rnberg! 
20071: llndebetouw fc Shulill2C I04a.b: ShuU & Slavi niTToM 
Slavin et alj 119931 iSlavin &: Coxi ,1992i: ,Borkowski et alj 
1990; Edg ar fc Chevalierill986D 

One of several possible production mechanisms for the 
high ions is turbulent mixing in places where hot and 
cool gas flow past each other. This idea is supported 
by the fact that the universe is also rich in hot. X-ray 
emitting gas and that the regions containing high ions 
are sometimes correlated with hotter, X-ray emitting re- 
gi ons. Early analytic models of mixing layers were made 
bv lBegelman fc Fabian! (1990). They estimated the tem- 
perature of the mixed layer as the geometrical density- 
weighted mean of the hot and cool gas temperatures un- 
der the assumption that two gases mix due to turbu- 
lence and that m ixing is efficient. iSlavin eFall (Il99l 
further developed Bcge lman fc FabianI (|1990l )'s idea and 
analytically calculated the emission spectra and column 
densities of high ions under the assumption that mixing 
reaches steady state. Their calculated ratios between the 
quantities of different ions are often used as diagn ostics 
for obser vations. The results oflSlavin et al.l ()1993l ) were 
tested bv lEsquivel et al.l ()2006[ ) who used 3-D magneto- 
hydrodynamic (MHD) simulations. They assumed coUi- 
sional ionization equilibrium (GIF) for their column den- 
sity calculations. They simulated the mixing layer on a 
10 pc wide computational domain for durations up to 3 
Myr, but they f ound that the i r simu lations did not reach 
steady state as ISlavin et al.l (|1993( ) had assumed. Al- 
though their calculated ion ratios can differ significantly 
with time and sightli ne, their ratios are more similar to 
the analytic results of lSlavin et al.l (|1993D than the aver- 
age ratios calculated for other phenomena such as radia- 
tive cooling. 

The next logical step in the progression is to calcu- 
late the ion content using non-equilibrium ionization al- 
gorithms. In this paper, we do that. We run detailed 
hydrodynamic simulations that include non-equilibrium 
ionization (NFI) calculations of interesting ions. We also 
perform CIE calculations for comparison. In order to 
offset the memory demands of the NEI calculations, we 
reduce the spatial dimensions from three to two. Our 
NEI results provide better predictions for high ion col- 
umn densities and ion ratios and can be used as updated 
diagnostics for comparison with observations. Our sim- 
ulations reveal that NEI calculations predict more high 
ions than CIE calculations, although the ion ratios do 
not change dramatically. We examine the reason for the 
difference between NEI and CIE by looking into the dis- 
tribution of ionization levels of interesting atoms fi J3.1.3p . 

We also consider various physical conditions that form 
turbulent mixing layers and find that the detailed con- 
figuration of the mixing layer, such as the perturbation 
amplitude of the initial boundary between the hot and 
cool gas, does not affect the ion ratios as long as mixing is 
efficient. However, the initial velocity difference between 
the hot and the cool gas and the temperature of the hot 
gas affect the ion ratios. When the initial speed has a 
slower value (50 km s~^) and the hot gas has a higher 



temperature (3 x 10^ K), the ion ratios show different be- 
havior than when the initial speed is 100 km s^^ and the 
hot gas is 1 X 10® K in temperature. These trends may 
be used as diagnostics that shed light on the conditions 
of the mixing gas fi )3.4.4[) . 

When c ompared with a nalytic turbulent mixing mod- 
els by Sl avin et al.l ()1993( ). our model simulations both 
from NEI and CIE calculations produce somewhat 
smaller ratios of C IV column density to O VI col- 
umn density, N(C IV)/N(0 VI), and/or somewhat larger 
N(N V)/N(0 VI) ratios. Our io n ratios are mo r e sim ilar 
to the CIE simulations done bv lEsauivel et al] ()2006D in 
which radiative cooling was allowed. We find that our 
average simulated ion ratios are close to the average val- 
ues from halo observations. Individual observations can 
vary greatly from the average. Some of the variation in 
the observed ratios may be due to line of sight geometry 
or age of the mixing layer, given that the ion ratios calcu- 
lated for individual sightlines in our domain vary greatly 
with time and sightline location. We also compare our 
simulated ion ratios with those observed in Complex G 
and find that our model simulations (and other turbulent 
mixing models) are more likely to produce the observed 
high ion ratios in Complex G when the low metallicity of 
Complex G is considered. 

This paper is organized as follows. The next section 
provides the numerical methods and physical parameters 
used in our simulations. Section 3 explains the results 
of our numerical study and Section 4 compares them 
with observations of the halo and Complex C. Section 
5 presents the summary. 

2. NUMERICAL METHODS AND PHYSICAL PARAMETERS 

We use FLASH version 2.5 for our simulations 
(iFrvxell et al. 2000). We include radiative cooling by us- 
ing relevant modules in FLASH. We also use the FLASH 
NEI module to track the degree of ionization of the car- 
bon, nitrogen, and oxygen atoms in the gas in each zone 
in a time dependent fashion. In FLASH, each timestep's 
NEI calculation is done in two pieces. In the first piece, 
which is done as a part of the hydrodynamic update, the 
mass density of each ion is updated via the mass conser- 
vation equation. During the hydrodynamic update, the 
total mass density and temperature are also updated. 
They will be used as inputs for the second piece of the 
NEI calculation. In the second piece, which occurs after 
the hydrodynamic update, the populations of the ioniza- 
tion levels within the atoms (of carbon, nitrogen, and 
oxygen) are updated by solving sets of ordinary differ- 
ential equations which include the ionization and recom- 
bination rates for each ion. In FLASH, the ionization 
and recombination rates include the effects of collisional 
ionization, auto-ionization, radiative recombination, and 
dielectronic recombination. 

In the FLASH NEI calculations, we set the abun- 
dances of carbon, nitrogen, and oxygen to be con- 
sistent with those used i n the Raymond and Smith 
code (Ra vmond fc SmithI (|1977l ): downloaded from 
HEASARC), which is used for the CIE calculations. In 
our CIE calculations, the Raymond and Smith code cal- 
culates the populations of the ionization levels in the 
atoms in each zone as a function of the zone temperature 
reported by the hydrodynamic simulations. The ioniza- 
tion fractions multiplied by the elemental abundance and 
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TABLE 1 

Models 



Domain 



Model 



Maximum 
Refinement 



Hot Gas 





X 


y 


Level 




T 


Vx 




(pc) 


(pc) 




(cm~^) 


(K) 


(km s 


A 


[0, 100] 


[-250, 50] 


6 


10-4 


10*^ 





B 


[0, 100] 


[-250, 50] 


7b 


10-4 


10** 





C 


[0, 10] 


[-25, 5] 


6 


10-4 


10^ 





D 


[0, 100] 


[-250, 50] 


6 


10-4 


lO** 





E 


[0, 100] 


[-250, 50] 


6 


10-4 


10^ 





F 


[0, 100] 


[-250, 50] 


6 


i X 10-4 


3 X 10** f 






Initial Interface 



y = f{x) 



Cool Gas 



riH 
(cm-^) 



T 

(K) 



(km s-i) 



S, = (2.5pc) sin(^) 
S, = (2.5pc) sin(^) 
y = (0.25pc) sin(^) 
S, = (2.5pc) sin(^) 
y = (5.0pc) sin(^)<= 
j/ = (2.5pc) sin(^) 



0.1 
0.1 
0.1 
0.1 
0.1 
0.1 



10'^ 
10^ 
10^ 
10^ 
10^ 
10^ 



100 
100 
100 
50 <^ 
100 
100 



Reference simulation 
^ Higher spatial resolution by one refinement level (factor 2) than model A 

The computational domain is 1/10 of model A 

Cool gas has half initial speed of model A 
° Amplitude of initial interface between hot and cool gas is twice that of model A 
' The temperature of the hot gas is three million Kelvin 



the volume density reported by the hydrodynamic code 
yield the volume densities of ions for each zone. The 
column densities are obtained by integrating the volume 
densities along sightlines. 

It is worth mentioning that our calculations are ap- 
proximate because we do not calculate the effects of NEI 
ionization levels on the radiative energy loss rate. Cal- 
culating the total radiative loss rate from the sum of the 
loss rates from individual NEI ions would require very 
large computing resources. Instead, in our simulations, 
we use the CIE cooling curve calculated from the CIE 
cooling rates of all of the relevant elements. Assuming 
that the plasma radiatively cools according to the CIE 
cooling curve at each time step of the simulation allows 
us to save significant computing resources. 

In order to verify the validity of this approximation, 
more complete future studies would be required, which 
compare simulations using NEI cooling with simulations 
using CIE cooling. The currently available comparison 
studies between CIE and NEI cooling rates show oppo- 
site trends depending upon whether the gas is in the pro- 
cess of ionizing (because its temperature has been raised) 
or is in the process of recombining (because its tem- 
perature has fallen! . [Sutherland k, Dopital ()1993 f) and 
ICnat fc SterlAeral (120071 ) calculated the NEI radiation 
rates from all relevant ions (as well as their ionization 
states) for cooling gas. Both studies showed that the 
NEI cooling rates are lower than the CIE cooling rates 
because the recombination of high stage ions in the cool- 
ing gas is delayed. In contrast, cool gas in the process of 
heating due to external heat sources shows delayed ion- 
ization and appears to have NEI cooling rates that are 
higher than the CIE cooling rate (jCnat et al.ll201Cih . In 
turbulent mixing layers, we can see both delayed recom- 
bination and delayed ionization in the mixed gas (this is 
verified by our simulations: see i jB.l.Sp . Thus the true 
cooling rate would be a complex combination of rates 
that sometimes exceed and sometimes fall below the CIE 
rate. 

We run our simulations in 2-D Cartesian coordinates 
for two reasons. Firs tly, the previous 3-D JVIHD study of 
lEsquivel et al.l (|2006[ ) showed that a 10 pc scale during 3 
Myr was not long enough to see efficient mixing of two 



gases. So, it is necessary to run larger scale simulations 
to a later time with significant spatial resolution. The 
mixing zone expands over time, so long duration runs 
must also have larger domains. Secondly, tracing ion 
fractions of carbon, nitrogen, and oxygen requires more 
computing resources. We find that 2-D hydrodynamic 
simulations with 100 pc x 300 pc domains running for 
several tens of IVIyr reveal the physical properties of high 
ions in mixing layers. In order to check the validity of 
our simulations, we also run a simulation with a 10 pc 
wide domain for 8 Myr. 

In this paper, we present the results of 6 simulations, 
labeled JVIodels A, B, C, D, E, and F. JModel A is our ref- 
erence simulation; each of the other models is made after 
varying one of IViodel A's parameters and runs in order to 
test the effect of that parameter. In IVIodel A, the com- 
putational domain is 100 pc x 300 pc and the maximum 
refinement level is 6. According to the FLASH adaptive 
mesh refinement (AJVIR) convention, this means that if 
our grid were fully refined, it would have 256 x 768 zones 
and the smallest zone would be (0.4 pc)^. At the begin- 
ning of the simulation, the hot gas {nn = 10"^ cm-'^, 
T = 10^ K) occupies the lower 5/6 of the domain and 
the cool gas {uh = 0.1 cm-'^, T = lO'^ K) occupies the 
remainder. The gas is initially in pressure balance with 
P/ks = 230 cm"^ K (considering that the number of he- 
lium atoms is 10% of that of hydrogen atoms). Note that 
this low thermal pressure is characteristic of the halo but 
the thermal pressure in the Galactic disk i s higher by at 



least an order of magnit ude (Ferrierc 199 jjJenkina[200l 
lCoxll2005l) . However, as lSlavin et 3111119931) pointed out. 



the column densities are independent of the thermal pres- 
sure although the emission intensity is proportional to 
the thermal pressure. The thermal pressure would also 
affect the speed at which the turbulent mixing layers de- 
velop. 

In order to create shear, the cool gas moves to the 
right (toward the positive x-direction) at 100 km s"^ 
throughout the simulation while the hot gas does not. 
Periodic boundary conditions on the left and right sides 
and outfiow boundary conditions on the top and bottom 
sides enable such motion. In order to seed the turbulence, 
we add curvature to the initial boundary between the 
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Fig. 1. — Model A: the top two panels depict the log of the hydrogen number density and the bottom two panels depict the log of the 
temperature. The leftmost plots in the first and third rows show the model at t = Myr. A time period of 5 Myr elapses between successive 
plots. 



cool gas and the hot gas. The shape of our boundary is 
y = (2.5pc) sin(^). 

Model B is a higher spatial resolution analog of Model 
A. It has a maximum refinement level of 7, which would 
correspond to 512 x 1536 zones if the domain were fully 
resolved. Model C has a smaller computational domain 
(1/10 in height and width), Model D has a smaller ini- 
tial speed difference (50 km s~^) between the cool gas 
and the hot gas, Model E has a larger ripple (maximum 
amplitude = 5 pc) in the boundary between the gases, 
and Model E has hotter (T = 3 x 10^ K) and less dense 
{uh = 5 X 10^** cm~'^) hot gas than Model A. Our model 



parameters are summarized in Table [TJ 

3. RESULTS 

3.1. Model A: Reference Simulation 

The sequences of snapshots of Model A presented in 
Figure [1] show the mixing zone between the hot and cool 
gas developing and growing in time. Five sequences are 
presented. They portray hydrogen number density, tem- 
perature, and indicators of the prevalence of C IV, N V, 
and O VI. The first snapshot in each sequence shows the 
model at i = Myr, when it contains only hot, rarefied 
gas in the lower portion of the grid and cool, denser gas 
in the upper portion of the grid. Between this epoch and 
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Fig. [T] continued. — Model A: population of C IV (top two panels) and N V (bottom two panels) in log scale as a function of time 
obtained from NEI calculations. In these and subsequent plots, the population of C IV, N V, or O VI refers to the fraction of ions in the 
specified ionization stage compared to the total number of ions of carbon, nitrogen, or oxygen, respectively. 

i = 5 Myr (shown as the next snapshot in each sequence), 
the relative motion between the hot and cool layers has 
stirred the gas, creating a zone of mixed, intermediate 
temperature gas between the hot and cool layers. Most 
of the mixing shown in our simulations occurs on the hot 
side of the hot-cool interface, where dense cool gas is en- 
trained and mixed into stationary hot, rarefied gas. The 
entrained gas expands due to the raising of its tempera- 
ture. Whatever hot gas becomes entrained into the cool 
gas compresses as its temperature falls and moves off of 



the grid at roughly 100 km s ~ 100 pc Myr 
appears only fleetingly in the simulations. 



and so 



When tendrils of cool gas first begin to intrude into 
and mix with the hot gas, the temperature of the en- 
trained gas rises immediately to the geometric mean. Its 
ionization state, however, does not equilibrate as quickly 
and thus the freshly entrained gas contains a mixture of 
poorly ionized gas that is in the process of ionizing and 
highly ionized gas that is in the process of recombining. 
As the next several snapshots (each of which is spaced 5 
Myr apart) show, the tendrils of intermediate tempera- 
ture mixed gas intrude into the hot gas as time goes by, 
increasing the depth of the mixed zone. 
The total depth of the mixed zone grows for the first 
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FlG.[T] CONTINUED. — Model A: population of O VI in log scale as a function of time obtained from NEI calculations. 



30 Myr then shrinks slightly, then begins to grow again 
around 50 Myr. Because the timescale of the growth 
spurts is a large fraction of the simulation time, it is not 
possible to conclude that the mixing layer has achieved 
a steady state, even after 80 Myr. Note that individual 
cool clouds in the Galaxy should be much shorter than 
8 kpc, the distance through which the hot and cool ma- 
terial in Model A slide past each other in 80 Myr. Thus, 
carrying out the simulations to larger times would not 
be realistic. 

Mixing plays an important role with regards to radia- 
tive cooling. Intermediate temperature (T ~ 10^ K) 
mixed gas cools far faster than hotter gas of the same 
density. Thus the mixing process seeds the hot gas with 
what it needs (cool gas) to cool quickly. Radiative cool- 
ing in the T = 1.5x10^ K, n = 6.8x10"^ cm"-^ mixed gas 
(such as that beneath the interface at i = 5 Myr) occurs 
on the timescale tcooi = | „a^) ~ (where A(T) is 

the cooling coefficient). As a result, this gas quickly cools 
to ^ lO"' K, the temperature at which cooling ceases to 
be applied by FLASH. As time progresses, neighboring 
mixed gas cools and joins it to make an ever-deepening 
layer of T ~ 10"* K gas. This layer contains, but is not es- 
pecially rich in C IV ions. It also contains N V and O VI 
ions, but in smaller portions to their atomic abundances. 
The hotter mixed gas is richer in all three of these ions. 

3.1.1. Calculation of Column Densities 

We calculate the average column densities of C IV, 
N V, and O VI ions along vertical sightlines, i.e., paral- 
lel to the y-axis through the Model A domain, and plot 
them at 1 Myr intervals in Figure [2] The plots present 
the column densities obtained using both the NEI and 
GIF algorithms, showing that the NFI algorithms yield 
much higher column densities than do the GIF algo- 
rithms. This will be examined in more detail in i )3.1.3l 
Galculating the column densities from the ion content of 



the entire domain is too GPU intensive, especially for 
GIF calculations. Instead, we find them by averaging 
the column densities of a large number of evenly spaced, 
infinitely narrow sightlines (256 sightlines for NFI and 
86 sightlines for GIF) spread across the domain. For 
the GIF column densities, we confirm that 86 sightlines 
are sufficient by comparing their average with that found 
from 256 sightlines for selected cases. Note that at any 
given time, the column density found along any given 
line of sight differs from those found along other lines of 
sight and thus differs from the average column density. 
The degree of variation is indicated, for example, in Fig- 
ure [U which shows column densities along each of 256 
sightlines through the grid at i = 30 Myr. Typical varia- 
tion ranges from as little as standard deviation/mean = 
16% for GIF O VI to as much as 52% for NFI G IV. The 
mean, median, standard deviation, and minimum and 
maximum values of the column density are summarized 
in Tabled] 

3.1.2. Evolution of Column Densities 

A common question asked while comparing simulations 
with observations or the results of analytic calculations 
is "which time frame should be used?" In order to un- 
derstand the answer, we discuss the time evolution of the 
high ion column densities here. 

The top-left panel of Figure [5] shows that the column 
densities of high stage ions are approximately at the 
beginning of the simulation when all of the gas in the do- 
main is either too cool or too hot to be rich in these ions. 
Once mixing begins, the high stage ion column densities, 
for both the NFI and GIF variants of Model A increase 
rapidly. Table |3] presents the amount by which the G IV, 
N V, and O VI column densities increase between t = 
and t = 10 Myr. Although the growth rates vary, the 
column densities of all three ions in both the NFI and 
GIF variants of Model A continue to grow until i w 20 to 
30 Myr, when the depth of the mixed zone peaks for the 
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Model A: Column Densities 



Model B: Column Densities 



Model C: Column Densities 




time (Myr) 

Model A: (NEI/CIE) Rotios 





time (Myr) 

Model B: (NEI/CIE) Rotios 



time (Myr) 

Model C: (NEI/CIE) Rotios 




VI ■ 
N V ■ 
C IV ■ 



Model A: N(C IV)/N(N V) 



time (Myr) 
Model B: N(C IV)/N(N V) 



VI ■ 
N V ■ 
C IV ■ 




time (Myr) 

Model C: N(C IV)/N(N V) 





20 40 60 80 

time (Myr) 



Model A: N(0 VI)/N(C IV) 




time (Myr) 
Model B: N(0 VI)/N(C IV) 



time (Myr) 
Model C: N(0 VI)/N{C IV) 






time (Myr) 
Model A: N(0 VI)/N{N V) 



time (Myr) 
Model B: N(0 VI)/N(N V) 



time (Myr) 
Model C: N(0 VI)/N(N V) 






time (Myr) time (Myr) time (Myr) 

Fig. 2. — Average high-stage ion column densities for all vertical sightlines in Model A (left column), Model B (middle column), and 
Model C (right column) as a function of time. Top panels: column densities of C IV (short dashed lines), N V (long dashed lines), and O VI 
(solid lines), where the dark lines are from NEI calculations and the gray lines are from CIE calculations. Second row of panels: ratios of 
NEI to CIE column densities for each ion. Third to fifth rows: ratios between the column densities of different ions; here the dark lines are 
from the NEI calculations and the gray lines are from the CIE calculations. 



first time. Between t ss 30 Myr and 80 Myr, the column 
densities of all three ions fluctuate. These fluctuations 
are also seen in the depth of the mixed zone in Figure [TJ 
However, the NEI C IV, in contrast with the NEI N V 
and O VI, climbs to a new plateau at i « 50 Myr. 

The average C IV column density found using NEI cal- 
culations continues to grow until late in the simulation 
because some of NEI C IV is in the somewhat dense, 
radiatively cooled region below the interface (see Figure 
[T]) and this region grows, although sporadically, during 
the simulation time. The gas density in the radiatively 
cooled region is larger than that in the warmer, actively 
mixing region below it. Thus, even though the fraction 
of carbon atoms in the C IV stage (in NEI simulations) 
is smaller in the cooled gas than in the warm gas (see 
Figure [T]), the number of carbon atoms in the cooled gas 
is significant. To demonstrate this point, we estimate 



the average column density of each high stage ion in the 
radiately cooled region and in the actively mixing region 
at i = 80 Myr, when the cooled region is most distinc- 
tive. We approximate the cooled region as that between 
y — —60 and y = pc and the actively mixing region 
as that between y = —200 and y — —60 pc. The col- 
umn densities are averaged over 256 sightlines for both 
the NEI and CIE calculations. The average column den- 
sities of each species in the cooled and actively mixing 
region are presented in Table S) It is shown that in NEI, 
the average C IV column density in the cooled region is 
47% as large as that in the mixing region and that the 
average volume density 1.94 x 10~^ cm""^) of C IV 
ions in the cooled region is even slightly higher than that 
(~ 1.79 X 10~^ cm~^) in the mixing region. 

The cooled region is poorer in N V and O VI than in 
C IV. Thus, the NEI column densities of N V and O VI in 



8 



TABLE 2 

Column Densities along Vertical Sightlines at t=30 Myr in Model A 



Ions 






NEI 








CIE 




mean 


median 
(10 


a 

12 cm-2) 


[min, max] 


mean 


median 
(10 


a 

12 cm-2) 


min, max] 


C IV 
N V 
O VI 


7.69 
1.49 
9.33 


7.36 
1.43 
9.29 


3.31 
0.35 
1.47 


[2.73, 16.17] 
[0.92, 2.41] 
[6.08, 12.19] 


1.56 
0.58 
4.18 


1.54 
0.54 
4.10 


0.81 
0.24 
0.80 


0.13, 3.87 
0.17, 1.32 
2.33, 6.98 





^ standard deviation 

Column Densities; Model A at t = 30 Myr 




Fig. 3. — Column densities of C IV (short dashed), N V (long 
dashed), and O VI (solid) ions along 256 vertical sightlines through 
Model A's domain at i = 30 Myr plotted as a function of position 
on the X-axis. Dark and gray lines are from NEI and CIE calcula- 
tions, respectively. 

the cooled region are not as significant as that of C IV. 
They are 19% and 17% of those in the mixing region, 
respectively. The C IV, N V, and O VI in the radiatively 
cooled zone principally derive from ions that were previ- 
ously hotter and are now in the process of recombining 
with free electrons. Since the CIE calculations ignore the 
plasma's history as they consider only the current tem- 
perature, they predict very few high stage ions in the 
radiatively cooled zone. In CIE, the column densities of 
all three ions in the cooled region are less than 5% of 
those in the mixing region. 

3.1.3. NEI vs. CIE 

As soon as mixing begins in our simulations, the NEI 
calculations predict greater numbers of high stage ions 
than the CIE calculations. (Because we assume CIE at 
the initial timestep, the NEI/CIE ratios at i = IVIyr 
are unity for all ions.) This is indicated by the second 
row of panels in Figure [21 and in particular, the left- 
most panel in that row, which plots the ratios of column 
densities predicted by the NEI calculations to the column 
densities predicted by the CIE calculations for IVIodel A. 

IMixing, which proceeds faster than ionization or re- 
combination, forces the gas out of equilibrium. For the 
first 10 Myr, most of the mixed gas is hotter than the 
CIE temperatures for C IV and N V, thus the CIE cal- 
culations predict low ion fractions for them. However, 
larger ion fractions and thus column densities of C IV 
and N V are found in the NEI calculations. These are 
due to previously-cool atoms that have lost a couple of 
electrons after mixing with hotter gas, but have not yet 
ionized to their equilibrium levels. Real observations of 
such gas should find wide absorption and emission fea- 
tures from this hot, 'underionized' gas. In addition, the 
profiles should also be widened by Doppler broadening 



in the turbulent velocity field. 

At their worst, the CIE calculations underpredict the 
C IV and N V column densities (relative to the NEI pre- 
dictions) by factors of ~ 8 and ^ 4, respectively, but 
this only occurs during the first ^ 12 IVIyr of the simula- 
tion. Subsequently, as the mixed zone develops a wider 
temperature profile, a radiatively cooled zone develops, 
and more ionization and recombination occurs, the CIE 
and NEI predictions begin to track each other with av- 
erage ratios of 4.6, 2.9, and 2.3 for C IV, N V, and O VI, 
respectively (see Table [5]). 

The effects of mixing, radiative cooling, ionization, and 
recombination can be seen in individual cells. Ideally, 
we could examine a single cell at various stages in its 
time evolution, from the moment after it first experi- 
enced mixing until long after it radiatively cooled. While 
this is not practical, it is possible to examine a vari- 
ety of cells at a single moment in time. Figure S] shows 
an expanded image of the domain and points out 6 cell 
locations. Cells 1-3 belong to the lower region of the 
computational domain, contain newly mixed gas on the 
boundary of the mixing zone, and are hotter. IVIost of the 
material in these cells came from the hot reservoir and 
minimal radiative cooling has occurred in them. Cells 4- 
5 are deeper in the interior of the mixing zone. Roughly 
half of their material came from the hot reservoir and 
half from the cool gas. These cells have radiated away 
significant fractions of their thermal energy and are now 
roughly 1/2 the temperature of the hot reservoir. Cell 6 
is nearest to the cool gas. Among all of our sample cells, 
cell 6 has mixed for the largest period of time, entrained 
the largest fraction of cool gas, and lost the largest frac- 
tion of its thermal energy to radiation. It has lost > 85% 
of its thermal energy to radiation and with a tempera- 
ture of ^ 15, 000 K, exemplifies the radiatively cooled 
zone at the 'base' of the mixed zone. 

Table [5] lists the locations (column 2) and hydrogen 
number densities (column 3) of these cells. The FLASH 
code allows us to trace the mass fraction of the initial 
cool and hot gas at each cell over time. The mass frac- 
tion of cell material that originally came from the hot 
reservoir {fhot, column 4) and the mass fraction from 
the cool reservoir {fcooi, column 5) are also listed in Ta- 
ble m From these mass fractions, we calculate the tem- 
perature that the gas would have had if radiative cooling 
had not occurred, T = That x fhot + Tcooi x fcooi, where 
That = 1.0 X 10^ K and T^ooi = 1.0 x 10^ K. This tempera- 
ture is tabulated in (column 6) in Table[6l The measured 
temperature at each cell from the simulation (T, column 
7) is the temperature that the mixed gas reaches after 
the radiative cooling is in effect. Comparing T (column 
6) and T (column 7) shows how significantly radiative 
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TABLE 3 

Column Density Increase at Early Times 



time period column density increase 

Model NEI/CIE [t^nin, tmax] 

(Myr) C IV N V O VI 



A 


NEI 


0, 10 


4.0 


0.56 


2.9 


A 


CIE 


0, 10 


0.49 


0.13 


1.1 


B 


NEI 


0, 10 


4.0 


0.59 


3.4 


B 


CIE 


0, 10 


0.76 


0.17 


1.2 


C 


NEI 


[0, 2] 


0.89 


0.19 


0.62 


C 


CIE 


[0, 2] 


0.17 


0.041 


0.38 


D 


NEI 


0, 10 


1.1 


0.18 


1.1 


E 


NEI 


0, 10 


4.4 


0.68 


4.4 


F 


NEI 


0, 10 


1.4 


0.2 


1.2 



TABLE 4 

Volume and Column Densities in Cooled and Mixing Layers of Model A 











C IV 


N V 


O VI 


NEI 


CIE 


NEI 


CIE 


NEI 


CIE 


Cooled ^ 


n ^ 


(cm" 




1.94 X 10-8 


2.69 X 10""^ 


1.22 X 10"9 


2.50 X 10-11 


4.36 X 10-9 


2.11 X 10-12 


layer 


N ^ 


(cm" 




3.60 X 10^2 


4.98 X 10^0 


2.26 X 10" 


4.63 X 10^ 


8.07 X 10" 


3.91 X 108 


Mixing 


n 


(cm~ 




1.79 X 10-8 


5.73 X 10-9 


2.80 X 10"'-' 


9.86 X 10-10 


1.07 X 10-8 


4.76 X 10-9 


layer 


N 


(cm" 




7.74 X 10^2 


2.48 X 10^2 


1.21 X 10^2 


4.26 X IQii 


4.62 X 10^2 


2.06 X 10i2 



y e [-60, 0] pc 
^ volume density of C IV, N V, and O VI 

column density obtained by multiplying volume density with the length of layer, 60 pc for cooled layer and 140 
pc for mixing layer 

y e [-200, -60] pc 



TABLE 5 

Ratios of Column Densities Calculated Using NEI to Those Calculated Using CIE 



C IV N V O VI 

Model 

mean median a [min, max] mean median cr ^ [min, max] mean median cr ^ [min, max] 



A 4.64 4.62 0.44 [3.82, 5.67] 2.91 2.88 0.34 [2.24, 3.82] 2.26 2.25 0.20 [1.83, 2.81] 

B ^ 4.30 4.27 0.54 [3.30, 5.53] 3.00 2.89 0.42 [2.28, 4.03] 2.29 2.32 0.25 [1.69, 2.83] 
C " 3.80 3.68 0.42 [3.21, 4.65] 2.73 2.66 0.34 [2.27, 3.50] 2.09 2.03 0.28 [1.73, 2.68] 



^ standard deviation 
^ averaged over t e [20,80] Myr 
averaged over t £ [6,8] Myr 

cooling has lowered the temperature of the mixed gas. 

As the temperatm'e drops (progressing from cell 1 to 
larger numbered cells) due to radiative cooling and ad- 
vection of cooler gas, the fraction of very high ions pre- 
dicted by CIE calculations decreases much faster than 
those predicted by NEI calculations. Because the NEI 
recombination rate lags the cooling rate so severely, we 
can find C IV in the 15, 000 K gas in cell 6, next to the 
cool gas. This is consistent with the sudden increase of 
column density of C IV at later times shown in the top- 
left panel of Figure H (123121) • 

Similarly, when cells advect cool gas, they gain low 
ions, which raise their low ion fractions. These ions are 
slow to ionize, causing the fractions of NEI C IV, N V, 
and O VI in cell 1 to be significantly greater than those 
predicted by CIE calculations. In cells 4 and 5, the NEI 
predictions for once, twice, and thrice ionized carbon, 
nitrogen, and oxygen are more similar to those of CIE 
calculations, but this is not solely due to recombinations. 
Instead, it occurs because radiative cooling has lowered 



the gas temperature enough for the CIE fractions of low 
ions to converge with the NEI ion fractions. Such gas 
is not actually in ionization equilibrium, as is apparent 
from the fact that the fractions of high ions greatly ex- 
ceed those predicted by CIE calculations. 

The hot reservoir is in or near collisional ionization 
equilibrium and thus mixing such gas into a cell raises 
the fraction of very highly ionized atoms in the cell. For 
this reason, cell 1 contains large numbers of hydrogen- 
like, helium-like, and fully stripped carbon, as indicated 
by the top-left panel of Figure [5l (Figure [5] shows the 
ion fractions, calculated from both the NEI and CIE al- 
gorithms, except for the NEI ion fractions of N VIII, 
O VIII, and O IX, which are not included in the plot al- 
though they were traced in our NEI calculations.) The 
greater prevalence of very high ions in the NEI case than 
in the CIE case is a sign of delayed recombination in 
NEI calculations; the recombination rate is slower than 
the mixing rate. 

The ionization behavior in our NEI simulations is 
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TABLE 6 

Cells for Ionization Level Calculations 





coordinates 


riH 


hot gas 


cool gas 


rp b 


T " 


cell 


(x,y) pc 


(10""' cm-3) 


mass fraction 


mass fraction 


(K) 


(K) 


1 


(24.9, -176.0) 


0.87 


0.76 


0.24 


760,000 


653,000 


2 


(34.3, -181.4) 


1.72 


0.65 


0.35 


650,000 


462,000 


3 


(40.7, -148.6) 


1.83 


0.60 


0.40 


600,000 


371,000 


4 


(59.5, -71.8) 


2.51 


0.44 


0.56 


441,000 


251,000 


5 


(78.0, -59.5) 


3.65 


0.36 


0.64 


361,000 


152,000 


6 


(61.9, -40.9) 


40.9 


0.11 


0.89 


110,890 


15,000 



^ Hydrogen number density 

^ Mixed temperature if there had been no radiative cooling. It is calculated from f = 
That X fhot + T^ool X fcool where T^ot = 1-0 x lO'' K, T^ool = 1-0 x 10^ K, and f^ot and 
fcool are mass fractions of hot and cool gas, respectively 
Temperature measured at the cell from the simulation 



consistent with pr evious studies with NEI calculations. 
I Ballet et al.l (|1986[ ) studied the evaporation of a spherical 
gas cloud and found that in NEI calculations, ionization 
to He-like stages is delayed. In their calculations, the gas 
is heated via conduction, while in our simulations, t he gas 
is heated due to mixing. iBoehringer fc HartauistI (I1987D 
added radiative cooling to the calculations o f IBallet et al.l 
(jl986[) and found a similar delayed ionization trend in the 
conductive interface of the evaporating cloud. Radiative 
cooling plays a more important role in the production of 
high ions when ho t (> a few times 10 ^ K) ga s cools radia- 
tively. Recently, iGnat fc Sternberg ()2007[) numerically 
calculated the ionization states and corresponding radia- 
tive cooling rates for the elements H, He, C, N, O, Ne, 
Mg, Si, S, and Fe. They included recent atomic data and 



Lag1 D Temperature (K) 
5a r 
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Cell 6 
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I 
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5.0 



4.5 



4.0 



3.5 



3.0 



25 50 75 1D0 
x{pc) 

Fig. 4. — Cell locations for ionization level calculations for car- 
bon, nitrogen, and oxygen. In this temperature map for Model A 
at t = 80 Myr, six cells are identified. The information for each cell 
is given in Table [6] and the ionization level calculations are showrn 
in Figure [5] 



investigated the effect of various metallicities. Their NEI 
calculations confirm that recombination is delayed rela- 
tive to the radiative cooling rate. Similar recombination 
lags have been sho wn in earlier NEI c a lculat ions, which 
are summarized in IGnat fc Sternberd (|2007f) . Previous 
CIE calculations in the context of radiatively cooled gas 
are also reviewed in IGnat &: Sternberg (|2007D . 

3.1.4. Ratios between Column Densities of Different Ions 

The third-left to fifth-left panels in Figure ^ 
show ratios between the column densities of differ- 
ent ions, N(C IV)/N(N V), N(0 VI)/N(C IV), and 
N(0 VI)/N(N V), respectively, calculated both from NEI 
and CIE for Model A. The ratios follow a similar trend 
as the column density evolution (top-left panel) and 
NEI/ CIE ratios for each ion (second-left panel) such that 
they start to stabilize around i = 20 Myr. Estimated 
mean, median, standard deviation, and minimum and 
maximum values of the column density ratios during the 
time period of t = 20 to t = 80 Myr for Model A are 
given in Table U\ 

The ratios of various ions' column densities are the 
most frequently used diagnostics for evaluating observa- 
tions because models that are based upon differing phys- 
ical processes predict different values. In the following 
sections, we will estimate the column density ratios be- 
tween different ions for all of our model simulations and 
summarize them in Table [71 Comparisons between our 
results and other models, including th e tu rbulent mix- 
ing calculations of lEsauivel et al.l ()2006D and lSlavin et alj 
([1993) will be given in g] 

3.2. Model B: Higher Resolution than Model A 

In Model B, the smallest cells are half the height and 
half the width of the smallest cells in Model A. All of 
the other parameters are the same as in Model A. The 
high ion column densities and ratios predicted for the 
Model B simulation are shown in the middle column of 
Figure [H The panels are placed in the same order as 
used for Model A (left column). A comparison between 
the left and middle column shows that the results of our 
high resolution simulation (Model B) are very similar to 
those of the reference simulation (Model A). The top- 
left and top-middle panels show that the evolution of 
column densities in Model B follows the same trends as 
Model A such that all column densities, both in NEI and 
CIE, continue to increase until f w 20 to 30 Myr and then 
fluctuate between t = 20 and t = 80 Myr. Similarly, the 
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Fig. 5. — Model A at t = 80 Myr: fraction of atoms in various ionization levels versus ionization level for three elements, carbon, nitrogen, 
and oxygen. The ionization levels are indicated by 1, 2, 3, etc on the abscissa, where, for example, 1 on a carbon plot refers to C I. Values 
from NEI calculations are indicated by triangles, while values from CIE calculations are indicated by squares. The NEI ion fractions of 
N VIII, O VIII, and O IX are not included in the plot although they were traced in our NEI calculations. Note that the CIE ion fractions 
at low temperatures are as small as the numerical precision in our CIE calculation code. The top row of panels pertains to cell 1 in Figure 
|4l while the subsequent panels pertain to cells 2 through 6 in Figure |4] Information about each zone is given in Table |6] 



column density ratios between NEI and CIE (second-left 
and second-middle panels) stabilize after t = 20 Myr for 
both models. As in Model A, the NEI calculations pre- 
dict more high ions than do the CIE calculations. Be- 
tween t — 20 and t — 80 Myr, the mean/standard devi- 
ation of the ratios of column densities calculated using 
NEI to the column densities calculated using CIE are 
4.30/0.54, 3.00/0.42, and 2.29/0.25 for C IV, N V, and 
O VI, respectively (Table [5]) . These values in Model B 
are close to those calculated in Model A ( §3.1.2|) . From 
the information in the third-middle to fifth-middle pan- 
els, we tabulate the mean, median, standard deviation, 
and minimum and maximum values of N(C IV) /N(N V), 
N(0 VI)/N(C IV), and N(0 VI)/N(N V), both for the 
NEI and CIE predictions. The values are given in Ta- 
ble [7] and confirm the similarity between Model B and 
Model A. 

Even though the overall results of Model B are very 
similar to those of Model A, we do find that resolution 
affects the temperature of the mixed gas and therefore 



the column densities at early times. Between t = and 
t = 10 Myr, the column densities of each ion in Model B 
(top-middle panel in Figure [2]) increase by 4.0 x lO^'^ 
(7.6 X 10"), 5.9 X 10" (1.7 x 10"), and 3.4 x lO^^ (1.2 x 
10^2) cm-2 foj. Q V, and O VI, respectively, in NEI 

(CIE) (see Table [HI). Ahhough t hese va lues for Model B 
are close to those from Model A f< j3.1.2p in both NEI and 
CIE, the increase in the column densities in Model B is 
slightly larger than that in Model A for all ions, which 
indicates that hot and cool gas mix slightly faster in the 
high resolution simulation between t = and t = 10 Myr. 
Because C IV in CIE traces the mixed gas around 10^ K, 
the gas around this temperature is most affected. It is 
produced much faster in the high resolution simulation. 

3.3. Model C: 1/10 Scale Simulation 

In Model C, we run the simulation for 1/lOth the time 
period (8 Myr) in a 1/lOth scale computational domain 
(10 X 30 pc) than in Model A. The purpose of perform- 
ing this model is to check whether the results of Model 
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NEI of Model D: Column Densities 



NEI of Model E: Colunnn Densities 



NEI of Model E: Colunnn Densities 





time (Myr) 
NEI of Model D; N(C IV)/N(N V) 



time (Myr) 
NEI of Model E: N{C IV)/N(N V) 





NEI of Model F: N{C IV)/N{N V) 

Model A 




time (Myr) 
of Model D: N(0 VI)/N(C IV) 



time (Myr) 
of Model E: N(0 VI)/N(C IV) 



time (Myr) 
NEI of Model F: N{0 VI)/N(C IV) 






time (Myr) 
of Model D: N(0 VI)/N{N V) 



time (Myr) 
of Model E; N(0 VI)/N(N V) 



time (Myr) 
NEI of Model F: N(0 VI)/N(N V) 




time (Myr) 





e (Myr) 



time (Myr) 



Fig. 6. — Column densities and ratios of column densities for Model D (left column), Model E (middle column), and Model F (right 
column) from NEI calculations compared with similar results for Model A. Top panels: column densities; gray lines are the Model A NEI 
results and dark lines are the NEI results of Models D, E, and F in the left, middle, and right panel, respectively. Second, third, and fourth 
panels: column density ratios of C IV/N V, O Vl/C IV, and O Vl/N V, respectively. 



A (our reference simulation) are still valid on a smaller 
scale and to study whether the difference between NEI 
and CIE begins at early times in the smaller domain. 
This is an important issue because we include radiative 
cooling which is not scale-invariant. If radiative cooling 
were not included, the smaller spatial scale simulation 
would be identical to the larger scale simulation because 
the spatial coordinates and time are proportional to each 
other in the hydrodynamics equations without any dis- 
sipative terms, making the calculations scale-invariant. 
For those researchers who are interested in the early evo- 
lution of the mixing zone, Model C can be seen as a 
close-up, with greater spatial resolution and smaller time 
periods. 

The column densities and ratios predicted for Model C 
are shown in the right column of Figure [21 in which the 
plots are placed in the same order (from top to bottom) 
as for Models A and B (left and middle column in Figure 
[2). At early times, the column densities of each ion in 
Model C (top-right panel) evolve in a similar manner 
as those in Models A and B, such that they continue 
to increase until i « 3 Myr, i.e., ^ 1/10 of the onset 
periods in Models A and B. The column density rises for 
each ion both in NEI and CIE for Model C are given in 
Table [31 These increases are close to those in Models A 



and B over 10 Myr, showing that the early evolution of 
the mixing layer does not depend on the scale as long as 
the radiative cooling rate is not faster than the mixing 
rate. 

There are ways in which Model C behaves more like 
Model A during its first 8 Myr than like a time-scaled 
version of Model A. In Model A, the depth of the mixed 
zone and the column densities of Li-like ions increase 
steadily for the first 8 Myr of the simulation. Simi- 
larly, they increase in Model C during this time period. 
But, if Model A and C acted like time-scaled version of 
each other, then the depth of the mixed zone and the 
column densities of Li-like ions would stall at 2 Myr 
(= 20/80 X 8 Myr) in Model C, which they do not. 
Model C's mixing zone grows so steadily that it over- 
flows the computational domain by 8 Myr. Model C, 
therefore, does not run long enough to develop a C IV- 
rich, radiatively cooled (T < 2 x lO'^ K) layer similar to 
that which appears in Model A around 10 Myr. As a 
result, most of the existing C IV in Model C resides in 
the mixed layer of actively cooling gas (T > 10^ K), the 
same region where O VI resides. The same can be said 
of Model A at early times. 
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TABLE 7 

Column Density Ratios between Different Ions 



Model 



N{C IV) / N(N V) 



N(0 VI) / N(C IV) 



N(0 VI) / N(N V) 



NEI / CIE 


mean 


median 


cr ^ 




min, max 




mean 


median 


a- 


min, max 


mean 


median 


a- 


min, max 


Ab 


NEI 


6.72 


6.75 


0.84 




5.04, 8.61 




0.70 


0.66 


0.21 


0.42, 1.22 


4.61 


4.38 


0.96 


3.14, 7.23 




CIE 


4.27 


4.07 


0.94 




2.39, 7.10 




1.46 


1.38 


0.48 


0.65, 2.96 


5.88 


5.72 


1.00 


4.26, 8.57 


B b 


NEI 


6.77 


7.05 


1.10 




4.20, 8.25 




0.68 


0.66 


0.21 


0.41, 1.23 


4.41 


4.30 


0.84 


3.18, 6.76 




CIE 


4.76 


4.73 


1.01 




3.01, 6.92 




1.28 


1.28 


0.43 


0.63, 2.32 


5.77 


5.71 


1.15 


3.68, 8.69 


C " 


NEI 


6.42 


6.43 


0.34 




5.90, 6.98 




0.59 


0.55 


0.08 


0.51, 0.72 


3.76 


3.69 


0.47 


3.09, 4.56 




CIE 


4.68 


4.57 


0.92 




3.49, 6.32 




1.08 


1.13 


0.19 


0.71, 1.31 


5.04 


5.10 


1.19 


3.25, 7.13 


Dd 


NEI 


8.17 


8.33 


0.94 


[6.08, 10.00] 


0.61 


0.49 


0.18 


0.42, 0.94 


4.81 


4.36 


1.03 


3.63, 6.85 


Eb 


NEI 


7.02 


7.17 


0.96 


[4.86, 8.69] 


0.67 


0.65 


0.20 


0.38, 1.25 


4.56 


4.33 


1.02 


3.17, 7.08 


F b 


NEI 


10.25 


10.37 


1.06 


[7.16, 12.55] 


0.55 


0.54 


0.11 


0.40, 0.96 


5.61 


5.56 


0.84 


4.39, 8.55 



^ standard deviation 
^ averaged over t G [20, 80] Myr 
^ averaged over t G [6, 8] Myr 
averaged over t S [30, 80] Myr 

3.4. Models D, E, and F: Modified Initial 
Configurations 

IModels D, E, and F are presented in order to address 
the following questions regarding the physical conditions 
of the mixing layer. Do different physical conditions af- 
fect the growth of mixing layers and do different physi- 
cal conditions affect the column density ratios such that 
the observed ratios can be used as diagnostics of mixing 
layers? Answering these questions, we find that modify- 
ing the initial amplitude of the disturbance between the 
hot and cool gases does not change the characteristics 
of mixing layer (see ij3A2] regarding IVIodel E). However, 
adjusting t he spee d difference between the hot and cool 
layers (see §3.4.11 regardi ng ]VIod el D) and the temper- 
ature of the hot gas (see §3.4.31 regarding IVIodel F) af- 
fect the column density ratios. Because NEI provides 
more realistic diagnostics than CIE, only the NEI results 
of these models are compared with those of IVIodel A, 
the reference simulation. Figure |6] shows the results of 
IVIodel D, E, and F in the left, middle, and right column, 
respectively, together with the NEI results of IModel A 
for comparison. 

3.4.1. Model D: Slow Speed 

In JVIodel D, the hot and cool gases slide past each 
other with a smaller speed (50 km s~^) than in IVIodel 
A (100 km s~^). IVIixing occurs more slowly in JVIodel D 
than in IVlodel A because mixing is caused by shear insta- 
bilities, which grow with the speed difference between the 
two bodies. As a result, the early-time high ion column 
densities build up more slowly (see the top-left panel of 
Figure [6] and Table |3]). Although slightly delayed, the 
column densities of Li-like ions in IVIodel D do approach 
those of Ivlodel A by the middle of the simulation time. 

The column densities from JVIodel D continue to in- 
crease until t « 30 JVIyr when they begin to fluctuate. 
We estimate the column density ratios between differ- 
ent ions for JVIodel D between t — 30 and t = 80 Myi 
and list them in Table [71 Compared with JVIodel A 
during a similar time frame, IVIodel D has a higher 
average N(C IV)/N(N V) ratio but a lower average 
N(0 VI)/N(C IV) ratio, while the N(0 VI)/N(N V) ra- 
tios are similar. The noticeable variations in the column 
density ratios due to only 50 km s~^ of velocity differ- 
ence suggest that the ratios may serve as diagnostics of 



velocity when the velocity difference between the hot and 
cool gas is large (see i i3.4.4p . 

3.4.2. Model E: Large Initial Amplitude 

The sine wave shape of the interface between the hot 
and cool gases is larger in JVIodel E than in JVIodel A, 
allowing us to test the effect of the interface's curvature. 
We find that JVIodel E's characteristics, including depth 
of mixing region, high ion column densities, and column 
density ratios, are similar to those of IVIodel A (see Figure 
[6] and Table [3]). These results show that the perturbation 
amplitudes in the initial interface between the hot and 
cool gas do not significantly affect the physical properties 
of the turbulent mixing layer as long as mixing is efficient. 

3.4.3. Model F: Hot Gas with Higtier Temperature 

In JVIodel F, the hot gas has a higher temperature (3 x 
10^ K) than in Model A. Because we require the thermal 
pressure of the initial hot layer to equal that of the cool 
gas layer, the hot gas in Model F has a smaller density 
(i X lO"'' H atoms cm"^) than that in Model A. 

The turbulent mixing layer in Model F evolves differ- 
ently than that in Model A; the high ion column densities 
increase at a slower pace and continue to rise throughout 
the simulation period (see Tableland Figure |6l). Fur- 
thermore the newly mixed gas is hotter and the radia- 
tively cooled gas (T < 2 x 10'* K) is shallower in Model F 
than in Model A because it takes longer for the hot- 
ter mixed gas in Model F to cool. We find that the 
high ion population is low in the hotter, newly mixed 
gas in Model F than in Model A. This effect seems to 
yield larger N(C IV)/N(N V) and N(0 VI)/N(N V) ratios 
than in Model A (see Figure [6] and Table [71). Model F's 
N(C IV)/N(N V) ratio is greater than that of Model A 
for nearly the entire simulation. Thus, this ratio, espe- 
cially if used in conjunction with the N(0 VI)/N(N V) 
ratio may serve as a diagnostic of the hot gas tempera- 
ture when interpreting observations of mixing gas. 

3.4.4. Diagnostics for Observations 

The column densities on sightlines that run perpendic- 
ular to the turbulent mixing layers are on the order of 
10" cm-2 for C IV, lO^^ cm-^ for N V, and 5 x lO^^ cm'^ 
for O VI in our model simulations (except Model C). In 
real observations, the sightlines probably pass through 
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the layers at sharper angles and intersect multiple mix- 
ing layers or multiple portions of curved mixing layers 
that surround individual clouds. Each of these effects 
increases the column densities by a multiplicative factor. 
When comparing with an observation, the model column 
density ratios are more useful than the model column 
densities because the multiplicative factor is unknown. 
In m.2\ we compare the ion ratios from our simulations 
with actual observations. Here, we discuss the use of 
column density ratios as diagnostics of the physical con- 
ditions producing the turbulent mixing layers. 

The results of Models D and F reveal that the ratios 
between the column densities of different ions can be used 
as diagnostics for the shear speed and hot gas temper- 
ature. If the N(C IV)/N(N V) and N(C IV)/N(0 VI) 
ratios along a sightline through a specific cloud are 
larger than those along sightlines through other clouds 
while the N(0 VI)/N(N V) ratio is similar along all 
of these sightlines, then it is likely that the mixing 
layer on the first cloud forms with a slower sheer speed 
(Model D). Another possible diagnostic comes from 
Model F. If observations along one sight line have higher 
N(C IV)/N(N V) and N(0 VI)/N(N V) ratios but lower 
N(0 VI)/N(C IV) ratios compared with observations of 
other clouds, it is possible to assume that the first mixing 
layer formed in hotter ambient gas than the others. 

Generally, more C IV is produced in the NEI calcula- 
tions of the models than in the CIE calculations, espe- 
cially in the radiatively cooled regions. This high popu- 
lation of C IV affects the column density ratios between 
different ions. However, C IV can also be produced by 
photoionization and is produced more easily than other 
high ions if there are nearby photoionizing sources be- 
cause the photoionization threshold of C IV is only 48 eV. 
Therefore, care must be taken when evaluating observa- 
tional results. 

4. DISCUSSION 

The column densities of high ions located in the Galac- 
tic halo are observable from the absorption lines in the 
spectra of halo stars and extragalactic objects. Here we 
compare observed high ion column densities and their 
ratios with those estimated from our model simulations. 
We also compare our results with predictions from other 
turbulent mixing models and from other sorts of models 
for high stage ions. However, there are some uncertain- 
ties in estimating the column densities from model calcu- 
lations. Therefore, we discuss these uncertainties before 
addressing how much they affect the results of model 
calculations and the comparisons with observations. 

4.1. Uncertainties in Model Calculations 

High ion column density predictions are strongly de- 
pendent on the assumed metallicity of the gas experienc- 
ing the modeled physical processes. However, the metal- 
licity in the Galactic halo is not well constrained and the 
uncertainty increases from the Milky Way to the inter- 
galactic medium and external galaxies. Therefore, the 
column density calculations vary according to the metal- 
licity used in the model. 

We find that most of previous model calculations used 
solar metallicity, which are sometimes better constrained 
by meteorite measurements than other methods, al- 
though the abundances of carbon, nitrogen, and oxy- 



TABLE 8 

Abundances Normalized to IAllenI (1197311 



Reference 


C " 


N 


O ^ 


C/N 


O/C 


O/N 


1 


1.34 


1.00 


1.12 


1.34 


0.83 


1.12 


2 


1.48 


1.07 


1.23 


1.38 


0.83 


1.15 


3 


1.10 


1.23 


1.29 


0.89 


1.17 


1.05 


4 


1.07 


1.02 


1.12 


1.05 


1.05 


1.10 


5 


0.74 


0.93 


0.74 


0.80 


1.00 


0.80 


6 


0.054 


0.013 


0.12 


4.15 


2.22 


9.23 


7 


0.74 


0.66 


0.69 


1.12 


0.93 


1.05 



References. — (1) Anders & Ebihara (1982) metal- 
licitie s use d in iBochrinKcr & Hartquist ( 1987) and 
in ISlavin £ Cox (1993) (2) Grcvessc (1984) used 
in jB orkowski ct al. (1990); Slavin ct al. ( 199 |) (3) 
lAnd crs & Grcvessc (1989) used in Sutherland & Dopital 
1199 3); Shclton (1998); Esquivcl ct al. ( 2006 !) (4) 
'Grcvessc ct al. (1996) used in Indebetouw & Shulj| 
(2004a). (5) Fox et al. (2004); adopted solar metallicity 
from HolwcKcr (2001); AUcnd e Prieto et al. (2002, 200y) 
(6) Fox et al. (2004); Complex C metallicity estimated 
along the PG 1259-1-593 sightline (7) .Asplund et all 
II2005I ); solar photospheric metallicity 
^ carbon, ni trogen, and oxygen abundance ratio to that in 
lAUenl 119731 ) for each reference. In AUe^ 119731) . carbon, 
nitrogen, and oxygen abundance is 3.31 X 10"*, 9.12 X 
10-5, and 6.61 X lO""* per hydrogen atom, respectively. 

gen still have large uncertainties, even in this method. 
Because metallicity measurements have been revised re- 
peatedly over the years, different model calculations used 
slightly different metallicities, especially for carbon, ni- 
trogen, and oxygen . For examp l e, the turbulent mix- 
ing calculations of ISlavin et al.l (I1993D 



- ■ used the solar 

photospheric metallicities of iGrevesse (IT981 and as- 
sumed that the cool gas has depleted metallicities (such 
that 50% of carbon, nitrogen, and oxygen are depleted 
into dust grains) while all dust grains are destroyed 
in the initial hot gas . In contrast, the CIE calcula- 
tions ([Beniamin et al.l 120011) in the numer i cal stu dy of 
turbulent mixing layers of lEsauivel et al.l ()2006D used 
the updated so lar photospheric metallicity data from 
lAnders fc Greve ssc (1989) which have differe nt abun- 
dances for carbon, nitrogen, and oxygen from iGrevess^ 
(fT984) (see Table [5] for the difference between these two 
metallicities) . 

For consistency between our NEI and CIE calcula - 
tions, we use the cosmic abundances of lAllenI ()1973f ) 
because they are the default abundances both in the 
FLASH NEI module and in the HEASARC version of 
the Raymond and Smith code, which we use for our 
CIE calculations. Note that the abundances of carbon, 
nitrogen, and oxygen from most of the reference s (be- 
fore th e year of 2000) are higher than those from lAllenl 
(119731) although recent rueasur e ments (.Asplun d et al.l 
[20051 I Allende Prieto et all [2001 [200lt [Holwegeiil200ll) 
measured lower abundances of these elements. Table [5] 
shows the ratios of solar metallicities from various ref- 
erences to the abundance from lAllenI ()1973() for carbon, 
nitrogen, and oxygen. The ratios for fractions between 
different atoms such as carbon/nitrogen, oxygen/carbon, 
and oxygen/nitrogen are also presented. They can be 
used to convert our column density ratios to those for dif- 
feren t metallicities. For e x ample , if the solar metallicities 
from lAnders fc Grevessd ()1989[ ) were desired instead of 
those from lAllenI (I1973D . then our column density ratio 
predictions would need to be multiplied by 0.89, 1.17, 
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and 1.05 for N(C IV)/N(N V), N(0 VI)/N(C IV), and 
N(0 VI)/N(N V), respectively. Note that this conversion 
for different abundances is approximate for the complete 
effects of abundances on the high ions (and their ratios) 
because abundances of metal ions also affect the cool- 
ing rates thus influencing the dynamical evolution of the 
gas containing the high ions. In optically thin gas, large 
abundances of metal ions increase the cooling rates be- 
cause the dominant cooling process is resonant line emis- 
sion from these metal ions. As mentioned in ij2l more 
complete future studies are required to address the full 
effects of abundances of all relevant ions on the cooling 
rates. 

Photoionization also significantly affects the column 
densities of high ions and the column density ratios be- 
tween different ions although we do not consider its effect 
in our simulations. Because the photoionization energy 
of C IV is much lower than those of N V and O VI, it 
is likely that photoionization preferentially enhance the 
C IV column densities. In this already mentioned 

in ij3.4.4[ it is not easy to distinguish observationally be- 
tween collisionally ionized C IV in the radiatively cooled 
part of the m i xing layer and photoionized C IV. Note that 
ISlavin et al.l (1199 3*) considered the effect of photoioniza- 
tion for the case where photons radiated from the mixed 
layer irradiate the cool gas. 

Besides the above two, there are two more uncertain- 
ties in dynamical models of turbulent mixing layers: time 
and sightline dependence. Column densities can be cal- 
culated along specific sightlines at specific times and can 
be directly compared with the observed ones. However, 
the computed column densities vary significantly with 
viewing angle and position, and in dynamical models 
also vary with time. For example. Figure [3] and Table 
[5] show column density variations between sightlines at a 
specific time {t = 30 Myr), while the top panels in Fig- 
ure [5] and ini show variations of sightline averaged column 
densities over time. In general, when these sources of 
variation are taken into account, the computed column 
densities from the dynamical models have wider ranges 
of predicted values than those from the static models. 

It is necessary to address again the issue of our simu- 
lated turbulent mixing layers reaching steady states, par- 
ticularly the question of whether the evolution of column 
densities and their ratios can be applied to interstellar 
and intergalactic clouds that are relatively small so travel 
the length of the cloud in less than 80 Myr (see ij3.ip . In 
our model simulations, the column densities and their 
ratios evolve until 80 Myr (except in Model C) which 
corresponds to a length of roughly 8 kpc for the 100 pc 
long mixing layers and a speed of 100 km s~^. This size 
may be too large for the realistic size of the clouds in 
the Milky Way although it is a possible size for some 
of HVCs relatively far away in the halo such as Com- 
plex C whose distance and area are measured as roughly 
10 kpc and 3 x 15 kpc^, respectively (jThom et al.ll2008D . 
In the case that the size of the cloud in the Milky Way 
is much smaller than 8 kpc, care needs to be taken when 
using the column densities of high ions in our simulations 
because the column densities continue to increase during 
the early times of the simulation. However, the ratios be- 
tween high ion column densities settle down quickly and 
do not vary significantly at later times implying that the 



ratios in Table [7] are applicable at early times as well and 
for small clouds. In addition, the Model C simulations 
which use a smaller domain and concentrate on the ear- 
lier times show that the ratios between high ion column 
densities settle down very quickly (around t ~ 2 Myr, 
which is similar to the cooling time scale in i j3.ip and re- 
main similar afterward, supporting the use of the values 
in Table [7] for earlier times (and for smaller clouds). 

Finally, we discuss the effects of some physical pro- 
cesses that we do not include in our simulations but 
could affect the results of our simulations. First, the 
results of our 2-D simulations for the NEI calcula- 
tions of high ion column densities could be different 
from more realistic 3-D simulations because it has been 
known that turbulent motion due to instabilitie s in 3- 
D is different from that in 2- D (jHussainI 119841: IBavlvl 
[19861: ICraik fc Cri minald [T986[) . The magnetic field is 
known to suppress turbulence both in 2-D and 3-D 
(iFrank etahl 119961: iJones et al.l [l997l: iJeong et "all 120001: 
IRvu et al.ll2Ci00ir These two effects may affect the col- 
umn densities of high ions because the region for the 
high ions' existence is where the turbulent motion oc- 
curs. However, as the comparison of our 2-D hydro- 
dy namic simulations w ith the 3-D MHD simulations 
of lEsauivel et al.l (|2006[ ) in the following section (g^j) 
shows, including the effects of the 3-D geometry and 
magnetic field would not significantly affect the ratios 
between high ion column densities. The turbulent mix- 
ing model would still distinguish itself from models with 
other physical processes regardless of these two effects 
(see Figure [T]). (Note that the CIE and NEI calcula- 
tions of our simulati ons predict similar ratios and that 
lEsauivel et al.l ()2006[ ) used only the CIE calculations for 
the high ion column densities.) 

The effects of thermal conduc tion for the tu rbulent 
mixing layers were discussed in lEsquivel et al.l ([2006). 
Following their estimation, our model parameters lead 
to a similar Spitzer thermal diffusion coefficient nsp ^ 
10^"* cm^ s~^ with Tmixed ~ 10^ K and Umixed 
10^^ cm^'^ but a larger turbulent diffusion coefficient 

Kturb ~ ^VturbLin] ~ 10^"^ Cm^ S"^ with L,nj ~ 100 pC 

and Vturb ^ 100 km s~^, where Linj is the energy injec- 
tion scale (approximated as the domain size) and vturt is 
the turbulent speed. In this estimation, we assume that 
the mixed gas has a typical temperature of T ^ 10^ K 
and a density of n ~ 10~^ cm""^ as i n lEsauivel et al.l 
(|2006[). As mentioned in lEsauivel et al.l (|2006D . the larger 
turbulent diffusion coefficient than the Spitzer thermal 
diffusion coefficient implies that the heat transfer is dom- 
inated by turbulence rather than thermal diffusion which 
has a smaller scale. However, the small scale diffusion 
process (including numerical diffusion) could affect our 
simulations, especially at early times. (The effect of spa- 
tial resolution that corresponds to the numerical diffusion 
is discussed in ^'3.2\) In our simulati ons as well as in the 
simulations of Esq uivel et al.l (j2006[) . the thermal diffu- 
sion coefficients at early times when the mixing zone is 
very shallow are larger than the above estimated values 
because the gas in the initial interface between hot and 
cool gas has a larger temperature gradient than the gas 
in the interfaces between mixed and hot gas or between 
mixed and cool gas at later times. Therefore, both in our 
and their hydrodynamic simulations, the heat transfer 
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Fig. 7.— Column density ra tios between different ions in log[N(N V)/N(0 VI)]-log[N(C IV)/N(0 VI)] s pace. Re-plotted from Fig. 1 
of llndebetouw &: Shulll l)2004al), with additions. Models: radiative cooling of Galactic fountain gas (friano/eg: IShapiro fc Beniam inI (1993); 
IBeniamin & Shapiro (1993)), turbule nt mixing layers ( asterisks; Slavin ct al. (1993)), conductive heating and evaporation of spherical and 
planar clouds (crosses; Bocliringcr & Hartauis9 |1987|) and Borkowski ct al. (1990), respectively), cooling supernova remnant shells (plus 
signs: rSlavin fcl3oxl ([1993); Shelton (1998)), and CIE (solid line; Sutherland & Dopita (1993)). The suite of average values from 3-D MHD 
turbulent mixing calculations of Esquivcl et al. ( 2006) are indicated by bright shaded regions (no magnetic field) and dark shaded regions 
(with magnetic field). Note that the shaded regions near the center are the results of their models with radiative cooling, while the shaded 
regions in the upper right are those without radiative cooling. The me dian values in Table [7] are plotted as filled (NEl) and empty (CIE) 
circles and the averaged value of halo observations along 34 si ghtli nes lllndebetouw SE Shullli2004h) is indicated by a filled star. The filled 
square is the result of the additional simulation mentioned in i|4.2l 



would be faster at e arly times with therm al conduction 
included. (Note that lEsquivel et all (j2006[) ran both hy- 
drodynamic and MHD simulations but did not include 
thermal conduction.) 

4.2. Comparison with Observations and Other Models 

The column density ratios between different ions for 
all of our model simulations are summarized in Table 
[71 Note that in our column density calculations, we use 
sightline geometries that are perpendicular to the ini- 
tial interface between hot and cool gas (i.e., perpendic- 
ular to the initial velocity vecto r of the cool gas) . This 
choice of sightl i nes is the same as lSlavin et al.l (jl993D and 
lEsquivel et al.l (|2006l ). Along these sightlines, the mea- 
sured velocities of each ion are so small that they can be 
directly compared with the halo observations (w w 0). 

The column density ratios in Table [7] are plotted in the 
log[N(N V)/N(0 VI)]-log[N(C IV)/N(0 VI)] space for 



the comparison with halo observations and other models 
(Figure ^ . The averaged value from the halo o bserva - 
tions along 34 sightlines in llndebetouw fc Shulll (|2004bf ) 
is indicated by a filled star in Figure [T] The median 
values in Table [7] are plotted as filled (NEI) and empty 
(CIE) circles. Error bars indicate minimum and maxi- 
mum values. 

The median values from our model simulations are 
with a fac tor of ~ 5 of the analytic turbulent mixing re- 
sults of Sl avin et al.l ( 19931). The 3-D hydr odynamic and 
MHD simulations of lEsquivel et al.l (j2006') in which ra- 
diative cooling was allowed produce similar column den- 
sity ratios as our simulations although they only used 
CIE calculations for the column density estimations (Fig- 
ure [7]) . Note that the bright and dark shaded regions 
near the center of Figure [7] represent their models that 
include radiative cooling, while the shaded regions in the 
upper right in Figure[7]represent the results of their mod- 
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els without any radiative cooling. (The bright and dark 
shaded regions indicate their models without magnetic 
field and with magnetic field, respectively.) The tur- 
bulent mixing phenomenon, evaluated by both analytic 
and numerical means yields larger N(N V)/N(0 VI) and 
N(C IV)/N(0 VI) ratios when compared with other phe- 
nomena such as radiative cooling, supernova remnants, 
and conductive heating (Figure [7]). 

The average ion ratios predicted by our models are 
similar to the average from halo observations, how- 
ever, our models produce more C IV and N V. From 
our model simulations, we find a trend that models 
with small velocity differences between the cool and 
hot gas (i.e. Model D) produce more C IV and we 
find that models with higher hot gas temperatures (i.e. 
Model F) produce more C IV and less N V. In or- 
der to confirm this trend and to find a case closer to 
the observed value, we run an additional NEI simula- 
tion with an initial speed of 150 km s~^ and a hot gas 
temperature of 2.0 x 10^ K. The result of this addi- 
tional simulation is plotted as a filled square in Fig- 
ure [7] with the same convention as the other models 
of our simulations (log[7V(N V)/7V(0 VI)] w -0.74 and 
log[A^(C IV)/A^(0 VI)] w 0.17). As expected, the filled 
square is shifted toward the star (observed value) from 
the data points of Mod el D and Model F. 

It was pointed out in Indebetouw fc Shulll (|2004bl ) that 
the collection of halo observations exhibits a wide range 
of column density ratios (their Fig. 4) that cannot be ex- 
plained by a single type of model. However, as discussed 
in i j4.11 the variations of column densities in dynamical 
models may explain the wide range of observed column 
density ratios. The error bars shown in Figure[7]are esti- 
mated only for the temporal variations after the column 
densities are averaged over all the sightlines at a given 
time. If the variations along sightlines are considered, the 
column density ratios estimated from our simulations are 
further scattered in the plot. 

Finally, it is interesting to compare the results of our 
simulations with the ratios of high ion column den- 
sities observed with high ve locities (i.e., highly ion- 

ized HVC s). For example, iCoUins et all (|2007h and 

iFox et ail (|2004[ ) measured the high ion ratios along 
three sightlines toward Complex C. The combined val- 
ues from their measurements are N(C IV)/N(0 VI) = 
0.44;° °^ and N(N V)/N(0 VI) = 0.19iEJ °? for Mrk 279, 
N(C IV)/N(0 VI) = 0.40t[!;!]^ and N(NV)/N(OVI) < 
0.11 for Mrk 876, and N(C IV)/N(0 VI) = 0.35t° g^ 
and N(NV)/N(OVI) < 0.07 for PG 1259-^593. (We 
choose the measured values when they are available from 
each reference article. When the measured values are 
different in different articles, we choose smaller values 
for the measured ion ratios and smaller upper limits for 
N(N V)/N(0 VI). Note that the difference between tow 
different values are not significant so that choosing larger 
values would not affect the following comparison.) These 
measurements would be close to the average halo obser- 
vations if they were plotted in Figure [T] However, the 
measured metallicities of carbon, nitrogen, and oxygen in 
Complex C are much lower than the solar metallicities. 
iFox et ail ()2004D measured the abundance of metals in 
Complex C along the sightline toward PG 1259-1-593 (see 
Table [8] for their measured carbon, nitrogen, and oxy- 



gen a bundances relative to lAllenI ()1973D ). iCollins et al.l 
(|2007D measured the average metallicity of Complex C 
as 0.13 solar metallicity based upon [O I/H I] measure- 
ments along 11 sightlines toward Complex C. They also 
found that there is much less nitrogen in Complex C by 
factor of 0.01 to 0.07 relative to solar abundance. 

When we apply different metallicities to the ion ratios 
calculated from our simulations, we need to shift the data 
points of our simulations according to Table |8l Apply- 
ing the measured abundances of metals in Complex C 
in Fox et al. (2004) to our calculations would shift all 
of our data points in Figure[7]by —0.35 and —0.97 along 
the N(C IV)/N(0 VI) and N(N V)/N(0 VI) axes, respec- 
tively. However, other data points from different models 
also need to be shifted by the same amount. This im- 
plies that our turbulent mixing model is more likely to 
explain the observations of highly ionized HVCs than are 
other physical phenomena. But, as mentioned before, 
the results presented in this paper pertain to sightlines 
that perpendicularly intersect the turbulent mixing lay- 
ers and so see nearly stationary velocities of calculated 
column densities. In a future study, we will investigate 
the velocity-resolved column densities, particularly fo- 
cusing on these HVC observations. 

5. SUMMARY 

We investigate the turbulent mixing layer model by 
running 2-D numerical simulations. Our simulations in- 
clude radiative cooling and NEI calculations. NEI cal- 
culations produce more high ions than CIE calculations: 
4.6, 2.9, and 2.3 times higher column densities for C IV, 
N V, and O VI, respectively for our standard model 
(Model A; Tables [D and [5]). We find that in NEI cal- 
culations, both ionization and recombination for these 
ions are delayed, resulting in more Li~like ions (C IV, 
N V, and O VI) than are present in the CIE calcula- 
tions (Figure [5] and §3.1.3|) . Many C IV ions, in partic- 
ular, are produced even in radiatively cooled mixed gas 
(T < 10"* K) at the base of the mixing zone because the 
recombination of C IV is slow in this relatively cool gas. 
These results are also valid in our other NEI simulations, 
namely those with higher resolution (Model B) and 1/10 
smaller computational domain (Model C). 

By changing the model parameters such as the shear 
speed, the initial amplitude of the interface, and the tem- 
perature of hot gas, we study the various configurations 
of the turbulent mixing layer. We find that more C IV is 
produced when the shear speed is smaller (Model D) and 
more C IV and O VI are produced when the hot reservoir 
has a higher temperature (Model F). The initial ampli- 
tude does not affect the column densities significantly 
(Model E). 

The study of various configurations shows that the col- 
umn densities and corresponding ratios between the col- 
umn density of different ions do not vary greatly when we 
vary the model parameters within the reasonable condi- 
tions for the mixing layers (a few times 10^ K for hot gas 
and a few hundreds km s~^ for the shear speed). When 
the results of our simulations are plotted together with 
previous model calculations in log[N(N V)/N(0 VI)]- 
log[N(C IV)/N(0 VI)] space (Figure El), they are consis - 
tent with tlie analytic estimations of Slavin et al.l (I1993D 
and the 3-D MHD calculations of lEsauivel et all (|2006D . 
All turbulent mixing layer calculations either analytic or 



18 



numerical predict more C IV and N V than the other 
models (radiative coohng, coohng SNR, and conductive 
heating) . 

We compare NEI and CIE predictions for a subset of 
our models, finding that the CIE calculations predict 
smaller values of N(C IV)/N(0 VI) and N(N V)/N(0 VI) 
(see the open circles for Models A, B, and C in the 
log[N(N V)/N(0 VI)]-log[N(C IV)/N(0 VI)] plot, Fig- 
ure [7]). However, the deviation of CIE from NEI is not 
as large as that between different phenomenological mod- 
els, especially between the turbulent mixing layer models 
and the other models, or as large as the variation between 
different sightlines or times in a single model. 

Because our NEI simulations are capable of calculat- 
ing the column densities of various ions along with ve- 
locity information, the velocity-resolved column densi- 



ties, especially in highly ionized HVCs, can be studied 
by using our simulations and comparing with obs erva- 
tions of previous studies such as lCoUins et al.l (|2007f ) and 
I Fox et ahl (j2004[ ) and future observations from STIS and 
COS newly installed on HST. We will investigate the 
velocity-resolved column densities in a future study. 

The FLASH code used in this work was in part de- 
veloped by the DOE-supported ASC/ Alliance Center for 
Astrophysical Thermonuclear Flashes at the University 
of Chicago. The simulations were performed at the Re- 
search Computing Center (RCC) of the University of 
Georgia. We appreciate the anonymous referee for his 
or her valuable comments on the radiative cooling and 
time scales. This work was supported through grant 
NNX09AD13G through the NASA ATPF program. 
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